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Abstract 

o 
o 

CN I We describe gluon-pair production from a space-time dependent chro- 

mofield via vacuum polarization within the framework of the background field 
method of QCD. The processes we consider are first order in the action. We 
derive the corresponding source terms for gluon pair production for situations 
that can be described by a 1+1 dimensional approach. Especially, we observe 

^ ' that within the range of applicability of our approach the principal contri- 
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bution to gluon production is included. Gluon production from a space-time 
dependent chromofield will play an important role in the production and evo- 

^^ . lution of the quark-gluon plasma in ultra relativistic heavy-ion collisions at 

§ : RHIC and LHC. 
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I. INTRODUCTION 

The quark-gluon plasma (QGP), a decon- 
fined state of matter, is predicted to exist at 
high temperatures around 200MeV and/or 
high energy densities around 2GeV / fm? [||. 
Apart from situations shortly after the big 
bang or inside a neutron star, heavy-ion col- 
liders could produce the necessary conditions 
for its existence 0. Even if this state of 
matter is then available in the laboratory, it 
will live for only such minute timespans and 
inside so small volumes that its direct de- 
tection is rendered impossible. This is the 
reason, why one is forced to look for in- 
direct signatures like J/'^ suppression [Q, 
strangeness enhancement |^,0, and electro- 
magnetic probes such as dileptons and pho- 
tons [§-|^. But, due to many uncertain- 
ties in the theoretical predictions and diffi- 
culties in the experiments, it is up to now 
still not possible to conclude about the exis- 
tence of the QGP based on observations of 
these indicators. That fact is mostly due 
to the deficient knowledge about the space- 
time evolution of the parton distribution in 
these reactions. Initially, they are in pro- 
nounced non-equilibrium. The way, the QGP 
relaxes towards equilibrium has got to be 
studied carefully. This can in principle be 
achieved through the solution of a relativis- 
tic non-abelian transport equation for quarks 
and gluons. Herefore, it is necessary to un- 
derstand the production of partons in ultra- 
relativistic heavy-ion collisions. There, the 
two colliding nuclei travel almost at the speed 
of light and are thus highly Lorentz con- 
tracted. When they pass through each other, 
a chromofield is formed between them due to 
the exchange of soft gluons |10|-[15| . This pic- 



tons collide with each other and proceede to- 
wards equilibrium. 

In all previous studies |]Tl],|14|,|T5[ , the pro- 
duction of the QGP is investigated by tak- 
ing parton production from a constant chro- 
mofield into account (which has been stud- 
ied in []TB|,|TP|). However, as seen in numeri- 
llJTH], the chromofield acquires a 



ture is the extension of the color flux-tube or 
string model widely used in high energy pp, 
e~^e~, and pA collisions |T6| , [T7[| . The chro- 
mofleld so formed polarizes the QCD vacuum 
and produces gg-pairs and gluons. These par- 



cal studies 

strong space-time dependence due to a com- 
bination of such effects as expansion, back- 
ground acceleration, color rotation, parton 
collision and parton production. In situations 
like these, parton production from a constant 
chromofield is not applicable and one has got 
to find the corresponding expression for a 
general space-time dependent field. 

e+e production from a space-time depen- 
dent chromofield is studied by Schwinger |]TB| . 
This is equivalent to computing the probabil- 
ity from the amplitude < /ci, fc2|'S'''^-'|0 > (de- 
noted by A^ -^ e+{ki)e-{k2)) with 5^^) ob- 
tained from the interaction lagrangian den- 
sity of a classical field and quantized Dirac 
fields 0]. Due to the same structure of the 
interaction lagrangian density, the produc- 
tion of a qq pair is similar to the e~^e~ case 
except for color factors |^. However, the 
computation of the probability for the pro- 
duction of gluons from a space-time depen- 
dent chromofield is greater than that for the 
production of fermionic, massive quarks and 
antiquarks. That is due to the larger phase 
space at the disposal for particles described 
by the adjoint representation. Hence, for a 
correct study of the production and equilibra- 
tion of the QGP, it is necessary to describe, 
how gluons are produced from a strongly vari- 
able chromofield. 

In [^ we have studied for the first time, 
the production of gluon pairs from a space- 
time dependent chromofield via vacuum po- 
larization. There, we derived corresponding 
source terms, discussed their range of appli- 
cability, and explicitely evaluated them for 
a time-dependent model field. However, the 



chromofield formed at RHIC and LHC is not 
only time dependent, but inhomogeneous in 
space and time. In the pre-equilibrium stage 
of the evolution of the QGP it is a good ap- 
proximation to consider a 1+1-dimensionally 
expanding system, i.e A depends on the time 
t and the longitudinal coordinate z. In this 
paper, we derive the source term for a t and z 
dependend chromofield. Also, for a compar- 
ison to Schwinger's formula for fermion-pair 
production (see Eq. 6.33 in fl^), we evalu- 
ate the total probability for gluon-pair pro- 
duction from a space-time dependent chro- 
mofield. In future, we hope to use the source 
terms derived in this paper for a t and z de- 
pendend chromofield in a non-abelian rela- 
tivistic transport equation to study the pro- 
duction and evolution of the QGP for situ- 
ations similar to those found at RHIC and 
LHC. 

Finally, we discuss the validity of our ap- 
proach to parton production, especially in 
view of an application to situations found 
at RHIC and LHC. As our approach is per- 
turbative, it is only capable of describing 
the production of particles with a momen- 
tum p larger than the product of the cou- 
pling constant g and the gauge field A. This 
is the principal contribution to the produc- 
tion of particles as long as gA is larger than 
the hadronic scale of QCD, i.e. Aqcd- So, 
we can use our method when the condition 
p > gA > Aqcd is satisfied. That means also 
that our approach breaks down as soon as the 
decaying field A reaches the value Aqcd/ 9 
from above. After this time tg^^^, a non- 
perturbative treatment of particle production 
is needed which is beyond the scope of this 
paper. In this paper, we investigate, whether 
the stage of particle production is already al- 
most over before we reach tend oi whether any 
significant contribution in the semi-hard sec- 
tor is still to be expected. 

The paper is organized as follows: In 
chapter II, we explain the setup of our 



method by deriving the probability for the 
production of gluon pairs to the lowest order 
in the effective action. In section HI, we give 
the source terms for a chromofield variable in 
1+1 dimensions. Chapter IV contains a brief 
summary. 

II. PROBABILITY FOR THE 
PRODUCTION OF GLUON PAIRS 

In this section, we will compute the prob- 
ability for the production of gluon pairs from 
a space-time dependent chromofield via vac- 
uum polarization. The process which con- 
tributes to the gluon pair production in lead- 
ing order of the action is diagramatically 
represented in Fig.(l). To evaluate these 
diagrams, we follow the background field 
method of QCD which, in a gauge invariant 
manner contains a classical background field 
and a quantum gluonic field simultaneously. 

In the background field method of QCD 



the gauge field A^ is split into two 
parts: ^4^-^71^ + Q^. Here, A^ is a classical 
background field which will not be quantized 
and Q^ is a quantum field representing the 
gluons. As the A-field is not quantized it does 
not occur in the functional integral and not 
in the coupling to the external source [dAi\ — >■ 
[dQ] and J A -^ JQ. 

In this method, a so called background 
field gauge: 

G« = d^Q^; + gr'^A'^Q'^ = {D^[A]Q^y (1) 

is chosen in a way as to make the total 
lagrangian C = Co + Cgf + ^fp with the 
lagrangian density for the gauge fields: 



Cg = —^F;^\A + Q]F'^^^\A + QI (2) 
with the contribution to the total la- 



grangian due to gauge fixing ||2^ 
Cgf ■ 



ya\2 



-,(G\ . 



(3) 



and the corresponding lagrangian density 
for the ghost fields 



for the three-vertex and: 



£ 



FP 



iD,[A]xi){D^[A + Q]x)a (4) 



M 



2A 



- I d%d^ki5^^\ki + k2-h-h 



gauge invariant with respect to type I 



gauge transformations |25l]. We evaluate the 
vacuum-polarization matrix element M =< 
^i^2|'S'^^''|0 > where S^^^ contains all inter- 
action terms of the Lagrangian density in- 
volving two Q-fields. For the diagram for the 
production of two gluons by coupling to the 
A-field once see Fig. (la). The Feynman rule 
for this vertex can be obtained from that part 
of the lagrangian density which involves one 
A field and two Q fields in the corresponding 
part of the lagrangian density and is given 
by: 



A^^{k,)A'^\k,)e''^{k,)e''P{k,)iy,A)tZ (9) 

for the four- vertex. Here A{K) is the 
Fourier transformation of the classical field 
A{x) given by: 



A{K) 



1 



(27r)^ 



d xe 



HK-x 



A{x) 



(10) 



Note that the above amplitudes include 
all the weight factors needed in order to re- 
trieve the correct lagrangian density. The 
probability is obtained from the amplitudes 
by the relation: 



{Via 



\abd 



w^Ej 



9r''[-29,pK-9,,iki 



^2)m + '^9puKp] 



spin 



d^k^ d^k2 

(2^0W 12^0^2^ 



MM* 



:ii) 



For the production of two gluons by 
coupling to the background field twice see 
Fig. (lb). The Feynman rule for this vertex 
can be extracted from that part of the la- 
grangian density which involves two A fields 
and two Q fields in the corresponding part of 
the lagrangian density and is given by: 



(5) 

To obtain the correct physical gluon po- 
larizations in the final state we use: 



spin 



spin 



(12) 



/ T A N abed 
\^2A)pu\p 



-^9V"'T''i9,x9up-9,p9uX 
+r'''r'%9,u9xp - 9,x9up 

+r''r'\9pu9xp - 



These Feynman rules coincide with those 
derived in [^. The amphtude for the pro- 
duction of gluon-pairs from the processes 
shown in Fig.(l) is given by: 



for the spin-sum of the gluons and after- 
wards deduct the corresponding ghost contri- 
butions which are shown in Fig. (2). For the 
integration over the phase space of the glu- 
ons there is an additional factor of 1/2 be- 
cause these particles (not the ghosts) in the 
9pp9ux)\- (d) final state are identical. Evaluating Eq.(|rT]) 
for the gluons, not yet including the ghost 
corrections, we find: 



9pv9xp) 
9pp9ux) 



W^ = WiAAA + WiA,2A + W2AAA 



W2A,2A 

(13) 



M = MiA + M2A 



(7) 



with: 



where 



(27r) 



W, 



M,A=^-^ j d^K5^'\K -k,-k2 



1A,1A 
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d^K^i 



A^'^iK)e'"^ik,)e'^ik2)iV,A)Z (§) 



pup 



X [{A^iK) ■ A*''{K))K^ - 

-{A'^iK) ■ K){A*''{K) ■ K)], (14) 



w 



IA,2A 



W' 



"iigac. 



2AAA 



d^K 



d% 



and: 



4 J " " (27r)2 
and: 



FP\bd 
2A ) 



w. 



otsg- 



2A,2A 



16 



[{{A'^^{ks)A'\K-k3))x 
x{A*^'^\k',)A*'^'^\K-k',)), (16) 



W 



- / d'^k'id'^ki X 

x5'^^\ki + k2-k:i-ki) X 
xA'^^{k,)A^\k,){yilr^^''. 

The probability in this case is: 

d^ki d^k2 



(23) 



FP 



FP\*bd 



(27r)32A;? {2T:f2kl 



(M^n (^ ) 



with: 



T, 



fiXfi'X' 



y / rab,cd _i_ rad,cb\/ ra'b,c'd _i_ ra'd,c'b\ 



Evaluating Eg . (|2^) yields: 



+ i2r 



9fitJ.'9xx' 



+ 

;i7) 



W^^^ = W^S^ + WZa 



(24) 



(25) 



with: 



with: 



wr 



-^ / rf^i^ X 



/ 



ab,cd 



abx rxcd 



fabxf 



:ii 



It has to be noted that the area of integra- 
tion for K is limited to (K^ > and K° > 0, 
because real gluon pairs are to be produced. 
Further, despite of a factor of i, Eq.(|^) is 
real. 

Now, we calculate the ghost part. The 
vertices involving two ghosts and one classical 
field and two ghosts and two classical fields 
respectively can be read from the lagrangian 
density and are given by: 



1A,1A 

X [{A''{K)-A*^{K))K^- 
-{A%K)-K){A*''{K)-K)] 



(26) 



and: 



W^^ 



as9 



2A,2A 



d'^Kx 



32(27r)4 
X {A\K) * A%K)){A'''{K) * A^'{K)) x 

/ _cab,cd _i_ rad,cb\/ ra'byc'd , ra'd,c'b\ {07\ 



{v,r)T = 9r\h - k,). 



and: 



/■j/-FP\abcd „'„2^ / fab,cd , fad,cb\ 

\V2A )fiX = -W 9f^X[J + J ) 



Note that W^^a = ^Hm = 0. To 
obtain the probability W for the produc- 
(19) tion of a real gg pair from a space-time de- 
pendent classical chromofield A, the result 
for the ghosts from Eq.([25|) has to be sub- 
tracted from the result for the gluons given 
in Eq.(ra). 



(20) 



The corresponding amplitudes for the 
ghosts are given by: 



(M 
with: 



FP\bd 



{M(^r + (M^^ir 



(21) 



{M[ir 



{2nf 



d^Kx 



x6^^\ki + k2 - K)A^^{K){VfX 



FP\abd 



f22) 



III. CHROMOFIELD IN 1-M 
DIMENSIONS 

We now give the source term for gluon 
production for a t and z dependent chro- 
mofield. This 1+1 dimensional approach is 
likely to propely describe the very early phase 
of a heavy ion-collision at RHIC and LHC 



The probability for particle production and 
the source term are connected by: 



99 



d^xd?k 64(27r) 



;A''^{xl)A''\xl) X 



W 



d^xd^k 



dW 



(2J 



X /rf2f^^«V'(^^)^c'A'^^^) ^ 



Note that of any expression given for a 
source term throughout the paper only the 
real part has got to be taken into account. 
The source term for the production of gluon 
pairs already corrected with the respective 
ghost contribution again can be composed 
analogously to Eq.(|l^). The general expres- 
sions for the gluon source terms given in |22| 
in 1+1 dimensions specialize to: 



e'^'^<'^-''^^^Ko{\kT\OfT;i 



/iA^'A'' 



(33) 



^n^'^^ 



Sg' 



d^xd^k 16(27r)4A;0 



A^^{xl) X 



X [8^^^/ - 3(/c^V + k' k',) - 



'kT 



X i^od^TlC) 



(29) 



where T differs from T by an additional 
factor of 2 in the second summand. The to- 
tal source term is the sum of Eq.(^), twice 
Eg. (|52[) , and Eg. (^3]) . Now, we ask the gues- 
tion, if the processes depicted in Figs. |1] and 
^ mainly only contribute to the particle mul- 
tiplicities inside the range of applicability of 
our approach. For this purpose we choose a 
special form of the field: 

A'^^ix) = Af^e-^'^'^eit + z)e{t - z)e{t), (34) 

which is confined to the inside of the light- 
cone. All other components are egual to zero. 
Many other forms of the field could have been 
chosen, but this choice is inspired by a nu- 
merical study presented in UTs). We define 
the time tf which gA needs to reach Aqcd'- 



where the transverse momentum compo- 
nents of k' are replaced by derivatives with 
respect to ixt- In principle, they can be cal- 
culated analytically, but the result is rather 
lengthy. Vectors with index L or T are always 
defined as follows: 

xl = (t, 0, 0, z) and fcy = (0, k^, ky, 0) (30) 

Additionally, we define the distance ( as: 



t 



f 



tnln 



gA 

A 



QCD 



(35) 



The remaining contributions are: 
3(^3 



(31) 



For to = 0.5/m, g = 1.5, Ai„ = 1.5GeV, 
and Aqcd = l50MeV tf is approximately 
egual to 1.3/m. The values for the de- 
cay time to, the coupling constant g, and 
the initial magnitude of the gauge field Ain 
are rough estimates for an ultra-relativistic 
heavy- ion collision at LHC (see [0). Now 
let us regard the accumulated density of pro- 
duced particles up to a time t: 



(mr2AlA 
99 



d^xd^k 8(27r)4fc0 



A 



aVVfi) 



/•* dW 
Jo d^xdt' 



(36) 



In our case this guantity is given by: 



2;z:l 



d^x' 



C 



^A''''{x^)A'^{x^)r'^''g^^> X 



2/12 



W 



g'A 



-2^ 



X e'^^-^^'^-^^'^K^{\kTmkL - tSL)x (32) 
and 



IGvrto 



[ll(e *o -e 'o) + 



-36/ALt^(e-^i-e-^^)]x 

xe{t + z)e{t-z)e{t). (37) 
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In Fig.d^) we have plotted the accumu- 
lated density w up to the time tf for the pa- 
rameters given above and for different values 
of the longitudinal coordinate z. One ob- 
serves right away that the principal contri- 
bution to the production of particles is from 
times significantly smaller than tf. To quan- 
tify this, let us concentrate on the accumu- 
lated number- density w at tf in the central 
region z = 0. We find: 



w 



IGvrto 



1^ 



^QCD^) + 



+ml{g'Al-kQCD% (3J 



From Eq.(^) one sees that for gAi^ » 
Aqcd the typical time to approach this end 
value of the density is of the order to- Which 
is by a factor of ln{gAin/AQCD) smaller than 
the time tf. So for the set of parameters 
chosen here, this factor is approximately 2.7. 
Again regarding Fig.@, we see that this ar- 
gumentation is well justified. From these 
facts also arises that the produced particle 
density is sharply peaked in the central region 
(see Fig.iQ)). This is because the size of the 
expanding system is still comparably small 
when most of the particles are produced. 



included within the range of applicability of 
our approach. Also, for a comparison with 
Schwinger's result for the total probability for 
the production of fermion anti-fermion pairs, 
we have evaluated the total probability for 
the production of real gluon pairs from an 
arbitrary space-time dependent chromofield. 
We plan to include the source terms derived 
in this paper in a non-abelian transport equa- 
tion with collisions between the partons taken 
into account. A self consistent solution of this 
equation would serve to study the production 
and equilibration of the quark-gluon plasma 
at RHIC and LHC. 
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IV. SUMMARY 



To summarize, we have derived the prob- 
ability for the production of a gluon pair 
from an arbitrary space-time dependent chro- 
mofield to the first oder in the effective action 
within the framework of the background field 
method of QCD. In order to make a connec- 
tion to the experimental situation, we have 
derived a source term for the production of 
gluons from a t and z dependent external 
chromofield. For this 1-1-1 dimensional de- 
scription of the plasma in the pre-equilibrium 
stage of an ultra-relativistic heavy-ion col- 
lision, we have observed that the principal 
contribution to the production of gluons is 
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Fig. 1 Feynman diagrams for the production of two gluons by coupling 
to the field A once or twice 
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Fig. 2 Feynman diagrams for the ghosts, corresponding to the gluon 
vertices. 
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Fig. 3 Accumulated density of particles w in GeV^ versus the elapsed 
time t in fm for different values of the longitudinal coordinate z. From 
top to bottom these are: z = 0.00 fm, z = 0.05 fm, z = 0.10 fm, z = 
0.15/m, z = 0.20/m, and z = 0.25/m. 
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Fig. 4 Accumulated density of particles w in GeV^ versus the longitu- 
dinal coordinate z in fm for different values of the elapsed time t in fm. 
From top to bottom these are: t = 1.3fm, t = 0.2/r7T,,and t = O.lfm. 
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